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Acoustic phonons can be coherently generated by ultrafast displacive screening of potential gradi¬ 
ents, often enhanced by the strong built-in piezoelectric fields, in wurtzite semiconductors. In such 
structures, transverse symmetry within the c plane hinders both the generation and detection of the 
transverse acoustic (TA) modes, and only longitudinal acoustic (LA) mode is generated. We show 
that even for c-GaN, the application of asymmetric potential distributions in the c plane can break 
the symmetry and selection rules, thus switching on the normally forbidden TA mode. This is in 
contrast to the LA mode, the strength of which varies with the symmetric potential distributions. 

By comparing transient differential reflectivity spectra in structures with and without asymmetric 
potential distributions, the role of the electrically attained anisotropy was further revealed by the 
digitized appearance of the TA mode, in clear contrast to the monotonic LA mode, and by modu¬ 
lations in the propagation velocities, optical birefringence, and geometrically varying sensitivities, 
the underlying mechanisms of which are modeled by electric-field-dependent perturbations of the 
dielectric tensors, incorporating the results of elastic modulations. 

BAGS numbers: 42.50.Wk, 63.20.-e, 78.20.hc, 78.35.+C 


I. INTRODUCTION 

Restructuring materials below the acoustic (AC) 
phonon mean free path or implementing variations in 
their elastic properties suggested new perspectives re¬ 
garding the controllability of the acoustic and ther¬ 
mal properties of crystals, such as thermal conductivity 
managemenfpii^ and charge transfer via AC pulses^, as 
well as the concept of phonon lasers^. Particularly in 
wurtzite crystals, where pseudomorphic strain caused by 
lattice mismatch at the interfaces leads to large internal 
polarization fields, photogenerated carriers screen out the 
strain-induced piezoelectric field, launching coherent AC 
motions^-—. This effect turned out to be very effective 
compared to thermal strain generation^ and deformation 
potential couplingi^. By virtue of strong electromechan¬ 
ical coupling, AC phonons in GaN-based piezoelectric 
heterostructures have been further investigated in terms 
of the strain dependenceii, tunability of zone-folded AC 
phononsi^, spatial modulation of AC wavesi^, and ter¬ 
ahertz (THz) electromagnetic radiation at acoustically 
mismatched interface s^^d^ . 

In this regard, GaN-based piezoelectric heterostruc¬ 
tures could serve as a versatile test bed for character¬ 
izing the AC functionalities with extraordinary genera¬ 
tion efficiency in the scheme of time-resolved ultrasonics. 
More specifically, (1) both zone-folded^ d^d^ and prop¬ 
agating AC phononsiiii^-— could be detected; (2) be¬ 
cause of the small absorption depth (<100 nm) compared 
to the conventional excitation spot size (>10 /im), the 
mode of the generated phonons depends on the growth 
direction; e.g., transverse acoustic (TA) and longitudi¬ 


nal acoustic (LA) phonons can appear simultaneously in 
the anisotropic planeiid^, or LA phonons alone can be 
generated in the isotropic c plan o^^d^d^ ; (3) both the 
deformation potential and displacive screening of electric 
fields, combined with the piezoelectricity, contribute to 
generating AC phonons having compositional amplitudes 
that depend on the crystal axis^. 

Owing to the well-established crystallinity of wurtzite 
semiconductors along the symmetry axis, where no shear 
mode can be excited by photoexcitation, and their poor 
detection sensitivity for TA phononsi^, previous studies 
were rather restricted to the LA properties among the 
three AC polarizations, although the TA modes are of 
great interest in the analysis of mechanical and optical 
properties and can provide nanoscale diagnostics with 
twice the temporal resolution because of their slow prop¬ 
agation velocities (approximately half that of the LA 
mode). To exploit the shear acoustic properties, sev¬ 
eral techniques based on broken lateral symmetry have 
been developed: tight focusing of pump beams within a 
spot size of 1^2 /im to reduce the AC directivity in the 
depth direction's, the transient interference patterns of 
two pumping beamsSi, LA-to-TA mode conversion at the 
isotropic-anisotropic crystal interfaceiS, and excitations 
on anisotropic crystal surfacesiSiSS. 

In addition, under an external reverse bias in GaN- 
based piezoelectric diodes, the amplitude of the LA 
mode propagating along the c axisiS and the local¬ 
ized TA mode along the non-c axis have been indepen¬ 
dently enhancedii. In most previous efforts, however, the 
growth direction was considered to be the crucial fac¬ 
tor determining the strength of the AC phonon modes 
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generated. From the viewpoint of prospective device ap¬ 
plications, it would be useful if an external bias could 
play a role in modifying the selection rules for generating 
and detecting AC phonons. We recently demonstrated 
such a switching scheme of TA phonons in a laterally bi¬ 
ased piezoelectric diodo^. The amplitude and switch¬ 
ing time of the TA mode could be electrically modu¬ 
lated. The structure with the electrically manipulated 
symmetry further revealed spectral shifts of the observed 
phonon frequencies due to modulation of the mechanical 
and optical constants. 

This work examines the relative roles of the electric 
potential distributions in switching and amplitude mod¬ 
ulation of AC phonons and presents the corresponding 
analytical expressions. For comparison, we measured the 
AC phonon properties of samples with or without an 
asymmetric potential distribution. The detailed theo¬ 
retical procedures are presented to describe AC modula¬ 
tion, including amplification and switching, on the basis 
of the electrically modulated elastic and dielectric proper¬ 
ties. The paper is organized as follows. We introduce the 
growth parameters, band structures, and strain-induced 
structural distortions of two contrasting samples, one un¬ 
der an asymmetrical potential distribution and the other 
under a symmetrical potential distribution, in Sec. IF 
Our experimental schemes are briefly described in Sec. 
III. In Sec. IV, we present the theoretical basis for elec¬ 
trical modulation of phononic properties, which is fur¬ 
ther classified into elastic tensor modulation for phonon 
generation and dielectric tensor modulation for optical 
detection of phonons. In Sec. V, we use our theoretical 
basis to explain the experimentally observed phononic 
spectra as a function of the external bias and dynam¬ 
ics. In Sec. VI, the role of crystal symmetry in modal 
detection is discussed in terms of the birefringence and 
detection sensitivity modulation under different poten¬ 
tial distributions. Finally, we summarize our findings 
and conclusions in Sec. VII. 


II. SAMPLES WITH DIFFERENT ELECTRIC 
POTENTIAL DISTRIBUTIONS 


Figure l(a,b) shows micrographs of two representative 
structures with distinct potential distributions. They 
have similar vertical structures with multiple quantum 
wells (QWs) placed within the intrinsic region of p-i-n 
diodes grown along the c axis, where the reverse applied 
bias compensates for the piezoelectric field^i; the funda¬ 
mental difference between them lies in the symmetry of 
the potential distributions. The sample under an asym¬ 
metric potential distribution (sample A hereafter), shown 
in Fig. 1(a), has a so-called interdigitated structure with 
narrow p and n electrodes that are laterally spaced by 
^100 pm along the x direction on top of a thin indium 
tin oxide (ITO) layer (^40 nm). In this way, an electric 
current / along the x axis was unevenly applied in sample 
A. In contrast, the contact electrodes of the sample under 




Sample S 




FIG. 1: Side and top views of (a) sample A and (b) sample 
S, where the blue dots along the x axis from the p contact 
to the n contact denote the laser excitation spots for the PL 
measurements. PL peak energy as a function of external bias 
in (c) sample A and (d) sample S. 


a symmetric potential distribution (sample S hereafter), 
shown in Fig. 1 (b), were much wider, with a thicker ITO 
layer (>150 nm) to evenly spread / in the x-y plane. 
The following layers were sequentially grown on sapphire 
substrates by metalorganic chemical vapor deposition in 
sample A (sample S): 3- (0.5-) /im-thick undoped GaN, 
2.5- (3-) /im-thick n-GaN, six (five) QW layers of 2 -nm- 
thick Ino.iGao.gN (1.5-nm-thick Ino. 25 Gao. 75 N) encased 
by seven (six) 8 - (7.5-) nm-thick GaN barriers, 120- (70- 
) nm-thick p-Alo. 05 Gao. 95 N (p-Alo. 15 Gao. 85 N), and 250- 
(140-) nm-thick p-GaN. The electron and hole concen¬ 
trations in sample A (sample S) were estimated to be 
about 7 X 10^^ (2 x 10^^) cm“^ and 8 x 10^^ cm“^, 
respectively. 

Figure l(c,d) shows the results of spatially re¬ 
solved photoluminescence (PL) measurements, where a 
frequency-doubled Ti:sapphire laser was used as an ex¬ 
citation source at 367.5 nm with a spot size of ^10 pm. 
The emission properties of the samples differed greatly. 
As the excitation spot moved away from the p electrodes, 
the PL peak energy was redshifted under the external 
bias Vext only in sample A. In the conventional case of 
sample S with an even current distribution, on the other 
hand, no notable energy shift was observed with x, as ex¬ 
pected. The PL energy shift is well understood in terms 
of the quantum-confined stark effecti^. A typical PL 
peak for spontaneous emission should be matched with 
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FIG. 2: Electric potential profiles of (a) sample A and (b) sample S near MQWs/n-GaN interface. Distributions of normal 
strain in (c) sample A and (d) sample S. Shear strain Szx distributions in (e) sample A and (f) sample S. Distortions in the 
lattice structures of (g) sample A and (h) sample S. 


the lowest-lying QW transition and determined on the 
basis of the net electric field in the QW, E'^. Hence, 
the redshift of the PL peak energy with x in sample A 
intuitively indicates that the actual magnitude of the ver¬ 
tically applied voltage decreased along x. To determine 
the energy band profiles in the x — z plane, E'^ was ex¬ 
tracted from the variational calculations of the PL peak 
energy. As E'^ and the net electric field in the barrier, 
E\^ are formulated using the built-in potential Vu^ the 
vertically applied voltage Va could be subsequently eval¬ 
uated as a function of x and z. The piezoelectric field 
Ep was set to 0.98 MV/cm in sample A and 3.7 MV/cm 
in sample S, as indicated by the solid lines in Fig. l(c,d). 

The resultant electric potential profiles, incorporating 
the exerted potential distributions of I 4 (x, z), are shown 
in Fig. 2 (a,b) for sample A at —4 V and sample S at 
—8 V, respectively. Following the procedures introduced 
above, the magnitude of E^^ (^0.8 MV/cm in sample A 
and ^0.5 MV/cm in sample S) reproduced the rapid 
potential variations along the 2 ; direction in both sam¬ 
ples. Further, the PL energy variations in Fig. 1(c) ulti¬ 
mately substantiated the increasing lateral electric field 
Ex = — in sample A, reaching ^0.66 kV/cm at —4 
V at the end of the depletion region. The magnitude 
oi Ex at —4 V was similar to that of E'^ in sample A 
but negligible in sample S. Accordingly, the tensor com¬ 
ponents of the mechanical strain are subject to distinct 
potential profiles in the piezoelectric structures exposed 
to the external bias. In particular, the additive lateral 
electric field Ex in sample A could induce shear strain, 
which is otherwise forbidden. 

According to the piezoelectric matrix, the strain com¬ 
ponents are connected by the electric field distributions: 


^xx 

^yy 

^zz 
^yz 
^zx 

where the piezoelectric strain constants are dsi = ds 2 = 
— 1.9 pmV“^, dss = 3.7 pmV“^, and dis = ^24 = —3.1 
pmV“^ according to the previously published data for 
wurtzite GaN with 6 mm symmetry^^. Then, the nor¬ 
mal (Czz) and shear strain components (Szx) in the x — z 
plane were calculated from the results in Fig. 1(c) [Fig. 
1 (d)] at —4 V for sample A (at —8 V for sample S) and 
are spatially mapped in Fig. 2(c-f), respectively. The 
signs of the normal strains {Szz) in Fig. 2(c,d) are op¬ 
posite to those of the in-plane strains {sxx and £yy)^ re¬ 
vealing positive Poisson ratios; the normal strains were 
compressive in the barriers and tensile in the QWs. On 
the other hand, the shear component (Szx) of sample A, 
shown in Fig. 2(e), became prominent throughout the 
biased region and was maximized near the end of the n- 
depletion region, in contrast to the suppressed values in 
sample S, shown in Fig. 2 (f), due to the absence of the 
lateral electric field. The strain-induced structural dis¬ 
tortions are illustrated schematically in Fig. 2(g,h). The 
volumetric deformation in Fig. 2 (h) for sample S, where 
the lattice constants are modified to c = co(l + Szz) and 
a = ao(l — i^Szz) from the strain-free values of cq and ao, 
implies that the structure is still laterally isotropic under 
an external bias Vext- Iri contrast, in sample A under a 
shear strain Szx^ monoclinic distortion that is geometri- 
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electricity. In conjunction with the types of structural 
deformations, the elastic and dielectric tensors are 
also classified into different symmetry groups. In this 
section, we develop mathematical formulations regarding 
the directional implications of the electric field com¬ 
ponents for modifying the acoustical properties of c-GaN. 


A. Electrical modulation of modal velocity 


FIG. 3: (a) Experimental scheme for time-resolved 

pump-probe and THz measurements using frequency-doubled 
Tiisapphire lasers, (b) Illustration of measurement geometry 
for probing AC modal dynamics. 

cally defined by an angle 0 to tilt the c axis^ breaks the 
optical and mechanical symmetry in the c plane in Fig. 
2(g). 


III. EXPERIMENTAL SCHEMES FOR 
MEASURING PHONON DYNAMICS 

Time-resolved pump-probe (PP) measurement and 
terahertz time-domain spectroscopy (THz-TDS) were 
conducted for various values of Vext oi* the angular pa¬ 
rameters (6>, (j)) to identify the role of the electric potential 
distributions in affecting both the detection and genera¬ 
tion of AC phonons. A pair of Ti-sapphire lasers, syn¬ 
chronized at 76 MHz and with a jitter (^113 fs) smaller 
than the pulse width (^250 fs), was used at either 367.5 
nm (for THz-TDS) or 367.5-375 nm (for PP measure¬ 
ment) in the reflective geometry, as shown in Fig. 3(a). 
For THz-TDS, each laser could be used as either a pump¬ 
ing source for THz generation or a signaling source (fixed 
at 800 nm) for detection via a Si lens coupled with a pho- 
toconductive antenna; the pump and probe excitation 
wavelengths were degenerate for PP measurement. The 
fluence of the pump beam for both PP measurement and 
THz-TDS was ^85 /iJ/cm^, and the ratio to that of the 
probe beam was 8:1. The laser spot size was maintained 
within ^20 /rm for the measurements. The angular pa¬ 
rameters of the measurements are further illustrated in 
Fig. 3(b). The incident angle 0 of the pump beam from 
the surface normal was fixed at 45° (for THz-TDS) or 
zero (for PP measurement), whereas the 0 value of the 
probe beam was set to less than 5° unless otherwise men¬ 
tioned (cf. Fig. 11-13). 


IV. EFFECTS OF ELECTRICALLY 
CONTROLLED CRYSTAL SYMMETRY ON AC 
PHONON DYNAMICS 

Applying an electric field parallel (perpendicular) 
to the 6-fold symmetry axis creates orthorhombic 
(monoclinic) distortion in a wurtzite structure via piezo¬ 


We first develop mathematical formulations explaining 
changes in the crystal symmetry depending on the poten¬ 
tial distributions and apply them to analytically express 
modal velocity variations. Under typical experimental 
conditions, the AC dynamics can be simply described by 
a loaded string model^i^^ assuming invariant crystal sym¬ 
metry, that is, fixed elastic coefficients. In our work, the 
influence of symmetry variation on the elastic properties 
is revealed to explain the experimental results of electri¬ 
cally controlled AC modal generations and modulations 
in the propagation velocities, which are experimentally 
estimated in Sec. V. 

As stated in Sec. H, the electric field components 
produce different types of strains. The effective elas¬ 
tic constants of a strained medium can be expressed as 
Cim = C^m + Cimn^n = 0^ + where is the 

constant for the undeformed medium, and Cimn is the 
third-order elastic (TOE) tensor that creates the pertur¬ 
bation ACim under the applied strain Sn. Such strain- 
induced elastic constant changes can alter the crystal 
symmetry's. For example, in c-GaN, the symmetry of 
ACim is still hexagonal for a normal strain Szz^ whereas 
the symmetry becomes orthorhombic for in-plane normal 
strains {sxx and Syy). When shear strains Szx and Szy are 
applied, the perturbed elastic tensors possess monoclinic 
symmetry with the diad axes parallel to the y and x direc¬ 
tions, respectively. In our case, the vertical and lateral 
electric fields induced normal {szz = dss^z) and shear 
strain {szx = di^Ex) via piezoelectricity. Therefore, the 
elastic tensor is expressed as 

O™ = CL + AC" + 

+ 0 ^ 3 ^ 33 ^. + ( 2 ) 

where AC^ (AC^) is the perturbed elastic tensor due 
to vertical (lateral) electric fields, representing hexago¬ 
nal (monoclinic) deformation, as illustrated in Fig. 2(h) 
[Fig. 2(g)]. The symmetry of the effective elastic ten¬ 
sor Cim consequently reveals that the structure could be 
monoclinically distorted^i^ under asymmetric potential 
distributions. AG wave propagation in such structures is 
described by the Christoffel equation, 

^klClrn^ran^n — PO ^^2 ’ (^) 

where \k is the atomic displacement velocity along the k 
axis, which determines the AG polarization. \/ki (V^n) 
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is the divergence (gradient) matrix operator, which is B. Electrically controlled modal detection 

given by 
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( 4 ) 


( 5 ) 


On the basis of the modified elastic constants, the ve¬ 
locities of AC modes propagating along the z axis with 
different polarizations are calculated to be 


vla 


sjch + C^sdssE. + 


\/p0 


(6) 


In this section, we mathematically formulate the re¬ 
lationship between the electrically manipulated crys¬ 
tal symmetry and optical sensitivities of AC phonons. 
An AC wave packet 77 ^( 2 :, t) propagating in a material 
causes a spatiotemporal change in the refractive index 
n by where i corresponds to the mode of 

AC wavepackets. In previous reports^ i^^i^^i^^ . Maxwells 
equations with phonon-induced refractive index modu¬ 
lations were developed for optically isotropic structures 
and used to obtain the differential reflectivity changes, 
AR/R^ in the probe pulse. Considering the dynamic 
Fabry-Perot (DFP) interference between the reflection 
from the surface, Vq = and the reflection from 

r]i{zA)^ Ari oc ^ 77 ^( 2 ;, t), AR/R is simply expressed as 
J dzrl^^rji(z,t)^ where is the detection sensitivity 
of r]i{zA) in an isotropic medium. In our experiments, 
however, the symmetry of the structures was externally 
controlled; hence, the detection mechanisms could differ 
from those previously reported. 


Jc2i + CgsdssE. - 

VTAll. - --, ( 7 ) 

and 


Therefore, we first characterize the refractive index of 
an electrically distorted GaN by and by taking into 
account the symmetry-dependent photoelasticity. Then, 
the propagation of the probe light along the 2 ; axis could 
be expressed by Maxwell’s equation: 


JcS, + C«,d„E, 

= V 7^ “v-7-^ ® 

Equations ( 6 - 8 ) represent that the modal velocities are 
externally manipulatable by applying electric fields. We 
also note that, without E^^ the TA modal velocity with 
its polarization parallel to the x axis {vta\\x) is the same 

as VTA\\y = ^ ^ 44 +^ 443 ^ 33 £^^ ^ Intcrcstiugly, for Ex (that 

is, for nonzero = di^Ex), Vta\\x < ^TA\\y, indicating 
that the lateral degeneracy of TA phonon velocities is 
lifted. 


-kl 

0 

o' 
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-kl 
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_ 0 _ 


where i^y^) is the electric field intensity of the probe 
light along the x (y) direction, the wave number k is 
<^pro 6 e/c, and e® is the unperturbed dielectric tensor of 
wurtzite GaN. The strain-induced dielectric tensor per¬ 
turbation, Ae, is described as 
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'pisdssEz pibdi^Ex 
PlsdssEz P25dl5Ex 
pssdssEz psbdi^Ex 

0 

psbdssEz + pAAdi^Ex 

0 


Pl^d^^Ez Plbdl^Ex 0 p35d33Ez PAAdi^Ex 

0 Pl3d33Ez -\-P2bdl3Ex 0 

P3bd33Ez + PA^dl^Ex 0 P2,2,d33Ez + P3bdl3Ex 


( 10 ) 


where the symmetry of the photoelastic tensor pim was 
determined by the structural distortion as pim = P^^n + 
^Pim = P^m Epimn^n based ou the unperturbed photoe¬ 


lastic constant p^^. The detailed form of pim is given in 
Appendix A. After solving the eigenvalue problem in Eq. 
(9), the refractive indices were calculated to be 


J 


Tighou — 


--[eo{ee -hP 33 d 33 Ez -\-p35d15Ex) + ee{pi 3 d 33 Ez -\-p15d15Ex) 


+ P33d33Ez + P35dl5Ex 

d33Ezdl5Ex{pi5P33 +^35^13 “ 2^35^44) + (^13^33 “ Pss) + (^15^35 - pIa)]^^^ (H) 

I- 


and 


nghu; = \/eo + PisdssEz + P25<^i5^cc, (12) 


where and eg are the ordinary and extraordinary com¬ 
ponents of e^, respectively. Similar to the case of AC 
waves, the refractive indices felt by the probe beam par¬ 
allel to either the x axis (righu;) or the y axis {fighoo) 
are the same without lateral electric fields; i.e., righoo 
=nghoj = ^/eo PisdssEz. In the particular case with 
only a lateral field but without vertical fields, nghoo 
\/^o+ PihdihEx + dl^E^{pi5ps5 - v\^pe and = 

E P 2 ^d\^Ex become nondegenerate, demonstrating 
the lateral anisotropy of the refractive index. 

The lateral anisotropy also affects the detection sensi¬ 
tivity Fi^'^ for normally incident probe light as 


r LA ^ 


dugh^ 

dps 


dTlghuj 


dps 



(13) 


M^ani 
r TA\\x 


OC 


dn, 


ghuj 


dr]5 





(pi 5 + 2disEx {pisPss 
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i^ani 
I- TA\\y 
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dn 


'Phoo 
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P55)Ae), 

(14) 

(15) 


where ps {ps) is the strain of the LA (x-polarized TA) 
mode. As a result, the amplitude of AR/R for the x- 
polarized TA mode depends on the probe polarization 


0 due to the distinct formulations between F^X\\x 
^TA||y’ where E^X\\j iiidicates F^X probe po¬ 

larization parallel to the j axis. 

The electrically attained optical anisotropy is further 
characterized in consideration of the distinguished de¬ 
pendence of Fi on the probe incidence angle 0. When 
AC modes propagate in an isotropic medium along the 
2 : axis, LA strain waves ps and TA strain waves 775 lo¬ 
cally modulate the dielectric tensor via the unperturbed 
photoelastic constants as 

'Pl3V3 

Pism 

A, ^ Pssm 

0 

pliPs 

. 0 . 

where the diagonal and off-diagonal components of Ae 
are separated by ps and 775 . Owing to the modally sep¬ 
arated tensor components, different modal sensitivities 
depending on 0 are expected as F^£X ^ 4/c'(l — 
and F£X ^ 2kxk'k{k' — k")~^{e^k'k" + k"‘^)~^p^^^ 
where k = ujprobelF ^x = ksm{0), k' = kcos{0), and 
/c" = e^k‘^ — k‘^. Conversely, the contribution of the 

diagonal (off-diagonal) components to AR/R can be rep¬ 
resented as F£X/Pi3 {E£X/p%a)' This explains why TA 
modes in an isotropic material cannot be detected with 
normal probe incidence, at which F£X becomes zero. 

In contrast, in an anisotropic medium, the dielectric 
tensor is modulated by combinations of ps and ps via the 


'pism 0 pliTis 
0 P? 3??3 0 

P44V5 0 


( 16 ) 































7 



FIG. 4: Illustration of the AR/R measurement scheme for 
two different AC wave packets, rjLA^zR) and rfTA^zR). 


perturbed photoelastic constants as 


'pism -\-pi5r]5 
PisVs -\-p25r}5 
P33V3 +P35^5 
0 

P35r}3 + P44P5 
0 


P13P3 + P15P5 0 P35P3 + ^44^75 

0 P13P3 P25P5 0 

P35r]3 + ^44^75 0 P33r]3 + ^35^5 


(17) 


Accordingly, measurements of the modal propagation are 
influenced by both diagonal and off-diagonal perturba- 
tions as F^a = AFLAlPi 3 )Pii + l3{Fff%lp%)p35 and 
Ft'a = AFla/Pi3)Pi 5 + l3{FTA/P4i)Pii- « (/?) is an em- 
pirical factor representing the mixed contribution of the 
diagonal (off-diagonal) perturbation in Ae to the signal 
amplitude of mode i via, e.g., the modified LA (TA) sen¬ 
sitivity due to the interaction between TA (LA) waves 
and probe beams. The sensitivities in the anisotropic re¬ 
gion are now simplified in terms of those in the isotropic 
region as 


value of a mode can influence the AT (Ar) value of other 
modes. 

Figure 4 shows the typical propagation of the strain 
pulses, the faster r]LA{z^t) and the slower r]TA{z^t)^ in¬ 
duced by the LA and TA mode in sample A, where the 
modal amplitudes are electrically controlled (cf. Fig. 5). 
The propagation of r]TA{z^ t), which lags behind rjLA{z, t) 
because its propagation velocity is ^50% slower, can be 
simply traced by measuring the DFP interference be¬ 
tween the reflection from the surface, ro, and the reflec¬ 
tion from r]TA{z,t), Avta—- However, the light reflec¬ 
tion from r]LA{z^t) is additionally modulated by a factor 
of (AFta)^, as r]TA{z,t) screens the probe beam inci¬ 
dent on and reflected from r]LA{z^t)^ which modifies the 
sensitivity function Fla and AR/R for the LA mode by 
the same amount. This acoustic screening effect on the 
modal sensitivity is experimentally verified in Sec. VI. 


V. ROLE OF SYMMETRY IN THE ACOUSTIC 
PHONON GENERATION 

In this section, we show that the application of an ex¬ 
ternal bias in the c plane can break the selection rules 
of isotropic c-GaN, switching on the normally forbidden 
TA mode. To explicitly show the differing roles of ver¬ 
tical and lateral electric fields in the generation process, 
a strain-dependent model of the AC wave equation was 
developed. Our experimental results further demonstrate 
that the modal propagation velocities are controllable un¬ 
der an external bias. In addition, the switching time of 
the TA mode can be effectively manipulated by changing 
the scale of the laterally biased region. 


A. Manipulation of AC mode selection rule 


r-ani\ 
r LA \ 
rani / 
i~ TAJ 


( aPl3/Pl3 I3P35/P4a\ (FTa\ 
\aPi5/Pi3 &PaTp\J \FTa) ’ 


(18) 


and are further compared with the experimental results 
in Sec. VI. 


1. Strain-based description of AC wave generation 

When coherent AC phonons are displacively initiated 
via the transient electric field screening by the photocar¬ 
riers, the relevant dynamic features can be described on 
the basis of a typical loaded-string equation^, 


C. Acoustic screening effect on modal sensitivity 

An AC wave packet locally perturbs the refractive in¬ 
dex of a material, inducing modulations of the light re¬ 
flection, Ar, and a corresponding transmission change, 
AF = — Ar, assuming that the phonon-induced absorp¬ 
tion change is negligible. In most ultrasonic studies with 
a single AC mode, these modulations, Ar and AF, have 
been investigated independently using different PP ge¬ 
ometries for reflectivity and transmission. However, for 
measurements of more than one AC mode, the com¬ 
bined effect of the modulations on the detection sensi¬ 
tivity should be taken into account, since the Ar (AF) 


d‘^Ui 2 

dt‘^ dz‘^ 


Si{x,z,t), 


(19) 


where i stands for each phonon mode (LA or TA), Ui 
represents the atomic displacement, and Vi is the modal 
velocity formulated in Sec. IV. The general form of driv¬ 
ing force Si{x, z, t) is expressed as 
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where Aa is the instantaneous stress change, eij is the 
piezoelectric stress constant, and AEj is the screened 
electric field intensity along the direction j. The driving 
force can also be represented in terms of the instanta¬ 
neous strain change Ae as 


Adxx 
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( 21 ) 


where is the unperturbed elastic tensor and [dij] 

is the piezoelectric strain tensor. Therefore, Aazz for 
Sla (=5..) is ( 2 |iC '?3 + where ^ -0.5 

represents the Poisson effect, and Adzx for St a {=Szx) 
is C^^Aszx- By assuming complete field screenings in the 
epicenters, which were confirmed by the saturated signal 
amplitudes in the fluence range shown in Fig. 6 (b), the 
modal driving force is 


Si{x,z,t) = —Ciei{x,z)H{t), ( 22 ) 

po 

where the roles of the externally applied strains and pho¬ 
tocarriers are clearly distinguished. The selection of the 
AC mode is determined by the strains, ela (^zz) and eta 
(^ zx), with the corresponding effective modal constants, 
Cla = + ^33 ^TA = C' 44 , respectively. The 

extra term 2^0^^ in Cla resulted from the contribu¬ 
tion of the in-plane strain changes {Aexx and AEyy) to 
the normal stress change Aazz- The instantaneous initi¬ 
ation of each mode by abrupt field screening is described 
by the Heaviside step function H(t) without regard to the 
mode index i. In this description, only the piezoelectric 
effect was considered because of its dominant role in the 
generation process in samples grown along the c axis^. 


2. Electrically controlled AC mode generations 

Ultrafast modification of the potential profiles by pho¬ 
tocarriers initiates the propagating AC waves; the probe 
laser waves encounter the modal AC waves with different 
positions and with spatial extents and are partially re¬ 
flected by or transmitted through the wave packet s^ 
Consequently, the multiple reflections of the probe laser 
waves form a DFP series, the oscillation frequency com¬ 
ponents of which are expressed as fi=2vincos{0t)/Xprobe^ 
where Ot is the angle of the probe transmission inside 
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FIG. 5: SWFFTs of the DFP interference from (a) sample 
A and (b) sample S under reverse bias from 0 to —5.5 V. 
The contour plots show spectral information (vertical axis) 
with respect to the time delay (horizontal axis). Integrated 
spectra over the entire temporal range for (c) sample A and 
(d) sample S. 


the materials and Xprobe is the wavelength of the probe 
beam. The transient oscillations in reflective PP mea¬ 
surements under Vext have accordingly been extracted 
from AR/R based on the sliding-window fast Fourier 
transform (SWFFT) spectra, as shown in Fig. 5(a,b). 
In sample A, not only /la at 119 GHz but also the novel 
spectral component at 70 GHz emerged under an exter¬ 
nal bias Vext- Ila matches well the vla value of ^7300 
m/s; this was further confirmed by the time of flight, 
which was revealed by an abrupt phase jump at the GaN- 
sapphire interfaceii, whereas the lower-frequency compo¬ 
nent (70 GHz) was ascribed to a new AC mode with vta 
(^4200 m/s^) and could not be explained by any previ¬ 
ous reports in isotropic structures. In clear contrast, sam¬ 
ple S without lateral fields showed a monochromatic spec¬ 
tral component the amplitude of which increased gradu¬ 
ally with the reverse external bias from 0 to —5.5 V. 

To verify that the 70 GHz component from sample A 
originated from the TA mode, the prominent frequencies 
measured with various values of Xprobe were converted 
into the corresponding velocities, as shown in Fig. 6 (a). 
The external bias was held at —4 V throughout the ex¬ 
periment, whereas the configurations of the pump and 
probe energies were degenerate except for one nondegen- 
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Sample A at -4 V 



Energy (eV) 


Excitation fluence (|jJ/cm) 


FIG. 6: (a) Modal velocities extracted from the DFP oscil¬ 

lations in sample A as a function of probe beam energy, (b) 
Spectral peak amplitudes of DFP oscillations as a function of 
pump fluence. 


erate case (3.32 and 1.66 eV, respectively) in Fig. 6(a). 
The high-frequency components (blue triangles) agreed 
well with the expected value of vla over the entire range. 
The 70 GHz component (red squares) also conformed to 
the Apro6e-dependent DFP characteristics of Jta at con¬ 
stant vta- 

In terms of the modal amplitudes, the bias-dependent 
changes in sample A were further explained by the differ¬ 
ent distributions of the normal and shear strains in Fig. 
2(c,e), respectively. To investigate the generation mecha¬ 
nism and relevant scale, the AC modal components in the 
Fourier-transformed spectra were traced as a function of 
the pump fluence at —4 V with a fixed laser wavelength 
of 367.5 nm, as shown in Fig. 6(b). The amplitudes for 
both AC modes were saturated around 80 /iJ/cm^, indi¬ 
cating that the generation mechanisms for both modes 
were dominated by transient field screening rather than 
deformation potential couplingi^. 

The strength of the screening electric field and the 
width of the epicenters can be estimated using a sim¬ 
ple virtual capacitor model^ formed of the accumulated 
photocarrier distributions. For coherent LA phonon 
generation, electric fields transiently induced by pho¬ 
tocarriers within multiple QWs (referred to as in-well 
screeningi^) play a major role. As Vext increases, how¬ 
ever, the reduced effective barrier widths result in con¬ 
siderable electronic tunneling, together with thermionic 
emission of holes over the lower valence band barriers^^. 
In this way, the charge densities of the serially con¬ 
nected virtual capacitors within the multiple QW region 
decrease; thus, the contribution of the QW carriers to 
AC mode generation also decreases. Therefore, the sup¬ 
pressed LA mode signal after —4 V in Fig. 5(a,c) implies 
that LA phonons originate from the QWs. 

At a large external bias Vext^ the macroscopic screen¬ 
ing, the spatial extent of which encompasses the in¬ 
trinsic and depletion regions (referred to as out-of-well 
screening^) may possibly become as important as the 
in-well screening. Under these conditions, the photo¬ 
carrier sweep-out under the vertical electric field (^0.8 
MV/cm at —4 V) leads to virtual capacitor formation. 


where the relevant sheet carrier densities are obtained by 
integrating the photocarrier distributions over the screen¬ 
ing distance. Using the photocarrier density observed un¬ 
der a pump fluence of 80 /iJ/cm^ and the high-frequency 
dielectric constant^, the required screening distance to 
compensate for the electric field of ^0.8 MV/cm was es¬ 
timated to be ^40 nm. Because of the significant sat¬ 
uration velocity difference between electrons and holes 
in GaN (a factor of the photoexcited elec¬ 

trons immediately accumulate at the boundary of the n- 
depletion region within ^0.4 ps, whereas the much slower 
holes take more than 5 ps to cross the entire macro¬ 
scopic region. Therefore, during the AC mode genera¬ 
tion period (^0.5 ps), a virtual capacitor scaled around 
40 could be spatially formed near the n-depletion 

region where shear strains are prominent, as shown in 
Fig. 2(e). In this regard, the TA signal, whose ampli¬ 
tude persistently grows with Vext despite QW tunneling, 
is thought to be generated from the n-depletion region. 


B. Symmetry-dependent AC mode dynamics 


1. Mode-dependent temporal lineshape analysis 

The propagation directions and spatial origins of the 
AC modes can be explicitly identified by tracing the am¬ 
plitudes and phase changes of the oscillatory components 
of AR/R. To explore different modal behaviors in sample 
A, the oscillatory signal (black line) at —4 V was band- 
pass-filtered at around the TA and LA modes (yellow 
scattered lines), as shown in Fig. 7(a). Because of the 
limited skin depth ^ of the probe beam (^700 nm), the 
envelope for a descending AC wave packet {r]f, where i 
indicates the phonon mode, LA or TA) would decay with 
time, whereas the envelope for an ascending wave packet 
{r]f) grows as r]f approaches the surfacei^. The ampli¬ 
tude change of each mode for 0 < t < in Fig. 7(a) could 
be fitted by considering both and rjf components (red 
and blue lines). 

When the values of are multiplied by the corre¬ 
sponding velocities, they are converted into the propa¬ 
gation distances of the AC phonons ascending from the 
epicenters to the surface (and vice versa for r]f). The 
distance for the LA mode (^415 nm) coincides well with 
the distance from the center of the intrinsic region to 
the surface (405 nm). Therefore, is considered to 
be generated mainly from the MQWs, as in conventional 
casesiii2^. In contrast, the TA mode traveled ^80 nm 
farther than the LA mode, which is approximately the 
sum of half the intrinsic region and n-depletion region 
widths, Thus, was considered to 

be launched mainly from the n-GaN depletion region, as 
can be further inferred from the shear strain distribution, 
which was maximized near the n-depletion region in Fig. 
2(e). For this reason, the slightly increased value of tta 
between —IV and —5.5 V (^ 20 ps) could be converted 
into the enlarged travel distance of 84 nm in accordance 
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Sample A at -4 V 




Sample A 


Wavelength (nm) 




Sample S 



Wavelength (nm) 


Wavelength (nm) 




FIG. 7: (a) Total oscillatory component in PP signal (black 

line) was filtered into the TA (52-88 GHz) and LA mode 
(108-144 GHz) frequencies. Inset illustrates the spatial ori¬ 
gins and propagations of the AG modes, (b) THz radiation 
as a function of time delay. Inset correlates the THz signal 
with the propagation of the ascending LA wave packet. 


with the expanded depletion region width at the surface 
and in the n-GaN region. This agreement between the 
transient digitized appearance and the travel path of the 
TA mode in Fig. 7(a) further confirmed that TA modes 
can be not only generated but also detected only in the 
laterally biased anisotropic region. This electrically ma¬ 
nipulated AC modal selection also agreed with the shear 
strain distributions in the biased region in Fig. 2(e). 

The initial deviations of the AC modes from the fit¬ 
ting curves are attributed to partial propagation of the 
surface-generated strain components toward the surface^ 
considering the relatively wide surface depletion region 
(SDR, ^150 nm)^. Without an ITO layer, the reflection 
of an AC mode from the air/p-GaN interface would be al¬ 
most 100%2iii^ without diffuse phonon scattering. In our 
case, however, the AC mode reflection from the ITO/p- 
GaN interface is estimated to be ^0.25 (0.19) for the 
LA (TA) mode, on the basis of the acoustic impedance 
mismatch between GaN^ and ITO^. 

The AC reflectivity at the p-GaN/ITO interface, r, to¬ 
gether with the spatial origins were further investigated 
by measuring the emitted THz waves during AC propaga¬ 
tions in Fig. 7(b). When an AC wave packet rji traverses 
a piezoelectric interface, transient polarization currents 
are induced as oc —{1/di^i — 

where di^i and d 2 ,i are the piezoelectric strain constants 
of the adjacent layers for mode i, as previously reported 
for LA wave s ^^4^ . Therefore, the transient THz emission 
patterns directly reflect the AC wave profiles crossing the 
interfaces at the position of Zint as illustrated in the in¬ 
set of Fig. 7(b); e.g., the temporal width of a THz pulse 
coincides with the spatial width of rji divided by Vi. 

To reduce complexity in the analysis, we focused on the 
propagation of the most intense strain component, 
dividing the THz signal into three temporal domains. In 


FIG. 8: Measurements of ellipsometric parameters (a,b) -0 
and (c,d) A from sample A and sample S for different Vext 
at wavelengths ranging from 350 to 380 nm. Vecct-induced 
modulations in dielectric constants for (e) sample A and (f) 
sample S. 

time domain / (0 < t < r^), travels from the mul¬ 
tiple QWs toward the ITO/p-GaN interface. The fast 
signal at the beginning is probably explained by surge 
photocurrent generation in the sample^. The signal at 
t ^ tla — 25 ps (denoted by arrow) is attributed to the 
interaction between 7^2 a p—GaN/AlGaN inter¬ 

face. In time domain II {ri < t < 2ri), the AC packet is 
reflected from the ITO/p-GaN interface with a reduced 
amplitude of r • r]^^. The amplitude ratio (^30%) of the 
signals at t ^ tla ± 25 ps (denoted by arrows) was in 
good agreement with the reflection of ^25% estimated 
from the AC impedance mismatch. The lingering slow 
oscillations in time domain III {2Ti < t < 3ri) are sus¬ 
pected to be generated by AC phonon-carrier interac¬ 
tions in the n-GaN sidei^. 


2. Electric-field-induced propagation velocity modulation 

Another intriguing feature of the electrically controlled 
symmetry is the spectral blueshifts of fi in both samples 
under increasing Vext in Fig. 5(c,d). To investigate the 
external-fleld-dependent dynamic properties of the AC 
modes using the information fi oc ViU, we first traced 
the refractive index changes with Vext by measuring the 
ellipsometric parameters, and A. These parameters are 
given by the ratio of the Fresnel reflection coefficients 
and Per for light polarized parallel and perpendicular to 
the incidence plane, respectively, as 

p = — = tan('0)e*^, (23) 

Per 

where the wavelength ranged from 350 to 380 nm, and 
the incidence angle was 75°. 

With the measured parameters and A in Fig. 8(a-d), 
the complex refractive indices h for both samples could 
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FIG. 9: FFT spectra for the ascending waves in (a) sample 
A and (b) sample S as a function of Vext- Spectral peak 
frequencies of AC modes in (c) sample A and (d) sample S 
under different Vext- Vecct-induced elastic constant variations 
in (e) sample A and (f) sample S. 


FIG. 10: AC mode amplitudes of (a) sample A and (b) 
sample S for 0 ranging from 0° to 90° in time domain I. 
Peak frequency for each AC mode in (c) sample A and (d) 
sample S as a function of 0. Refractive index ellipsoids for 
perturbed (at —5.5 V) and unperturbed (at 0 V) cases in (e) 
sample A and (f) sample S. 


be evaluated as 

„ _ y^l — 4sin^(^o) tan(i/;)e*^ + 2 tan(i/;)e*^ +tan^('0)e*^ 

^ sin(^o)“^ cos(^o)[l + tan^('0)e*^] ’ 

(24) 

where the real part of h is the refractive index n. The 
unperturbed dielectric constant 6*° of sample A (sample 
S) was evaluated to be 8.64 (9.12) at 367.5 nm; using 
this value, the modulation of the permittivity Ae{Vext) 
w as subsequent ly extracted based on the relation n = 
+ Ae{Vext)^ as shown in Fig. 8(e) [Fig. 8(f)] for 
sample A (sample S). 

On the other hand, the modal spectra of the DFP oscil¬ 
lations were traced without phase-change-induced errors 
by selecting temporal ranges for the FFT between Tij2 
and Ti, as shown in Fig. 9(a,b) for samples A and S, re¬ 
spectively. The peak frequencies of both the LA and TA 
modes in sample A increased by ^1.7% with increasing 
Vext ill Fig. 9(c). In contrast, the spectral shift of the 
LA mode in sample S in Fig. 9(d) was relatively negli¬ 
gible, ^0.25%. Then, inserting the measured n into fi 

under Vext', we extracted Vi equal to y/ C+ACffVext) ^ 

corresponding bias-dependent changes, AeiVext) in Fig. 
8(e,f) and ACiiVext) in Fig. 9(e,f), were further com¬ 
pared with the analytical expressions as Ae 

lipis +Pi5)di5E^, ACla - 

^ ^33 '^44 

and ACta — C^^^d^zEz - d^^El. Here, we in- 

'-^33 ^44 

corporated the optical birefringence and strain-induced 
nonlinear elasticity in the anisotropic region based on 
the photoelastic tensor component pij and the TOE co¬ 
efficient for the normal (shear) strain, "Fhe 

relatively moderate behavior of Jla with Vext in sample 
S implies that the shear strain contribution is dominant 


even in modifying vla^ as demonstrated in Eq. (6). 

VI. ROLE OF SYMMETRY IN THE 
DETECTION OF ACOUSTIC PHONONS 

In an optically isotropic medium such as c-GaN, di¬ 
electric tensor modulation due to shear strain can be de¬ 
tected only under oblique probe light incidenco^ii^. The 
shear strain in a laterally biased anisotropic medium, on 
the other hand, induces perturbation not only in the off- 
diagonal components, but also in the diagonal component 
of the dielectric tensor. Accordingly, the TA mode was 
optically detectable even when the incidence angle was 
normal to the surface, as represented by the digitized ap¬ 
pearance in Figs. 5(a) and 7. In this section, we experi¬ 
mentally verify the optical anisotropy and corresponding 
changes in the detection sensitivity by measuring the AC 
modal properties using different combinations of probe 
polarizations (j) and incidence angles 0. 

A. Electric-field-induced optical birefringence 

Regarding the electrically manipulated crystal sym¬ 
metry, Ex breaks the hexagonal symmetry of sample A 
via the TOE tensors into monoclinic symmetr}^. Con¬ 
cretely, the tilting of the principal axis away from the 
c axis by Ea,-induced monoclinic deformation could be 
illustrated by angular perturbations of the elastic and 
dielectric tensors in the x — y plane. 

Figure 10(a,b) shows the spectral amplitudes of DFP 
oscillations from both samples at a Vext value of —5.5 
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V, integrated over ri/2 < t < r^, as we rotate the probe 
polarization 0 from the x axis at ^ = 0. In sample A, the 
TA amplitude monotonically decreased by about 

50% with increasing 0, in clear contrast to the constant 
LA mode (A£^) in Fig. 10(a), indicating that the TA 
mode was partially polarized along the same direction 
as Ex. Figure 10(b) shows that in sample S, was 
invariant with respect to the probe polarization 0, as can 
be inferred from Eq. (13). 

Despite the large value of Vext^ the modal frequencies 
/la and /ta in sample A exhibited very similar increases 
with 0, as shown in Fig. 10(c), owing to the birefringence 
in the anisotropic region. In contrast, no notable change 
was observed for sample S, as shown in Fig. 10(d). As 
a function of 0, in this regard, the n value of sample 
A was extracted from either /la (blue squares) or /ta 
(red triangles) in Fig. 10(e). The dashed inner circle in 
Fig. 10(e) displays the isotropic refractive index ellipsoid 
without Ex, whereas the solid outer circle exhibits the 
values calculated using the F^a,-induced birefringence. In 
sample S, which has a symmetric potential distribution 
even under an applied Vext^ the Eea^t-induced birefrin¬ 
gence and nonlinear elasticity was not observed in Fig. 
10(f). For this reason, the hexagonal symmetry of the 
lattice structure should be preserved with righoo = righoj. 
The comparison of the samples therefore confirms the 
dominant role of the lateral potential gradient in manip¬ 
ulating the crystalline symmetry, concurrently with the 
modified elastic and optical coefficients. 


B. Optical characterization of anisotropy in sample 

A 

To further investigate the optical anisotropy in sam¬ 
ple A, we varied the probe polarization 0 at a fixed 0, 
now with an 5-polarized analyzer in front of the probe 
detector. When the probe beam is incident on the sam¬ 
ple at nonzero 0, the DFP interference for the TA mode 
is expected to be maximized at a 0 value of 45° in such 
a geometr}^. Since the PP signals became attenuated 
as they propagated downward toward the n-GaN side, 
the laser wavelength was further increased to 375 nm to 
ensure a deeper probe beam penetration depth ^ (^1.5 
/im). The ^-dependent DFP oscillations under —4 V are 
shown in Fig. ll(a,b) for normal and 55° probe inci¬ 
dence, respectively. To extract the spectral amplitudes 
of AC modes propagating in regions of different symme¬ 
try, we performed FFTs in two distinct temporal ranges, 
that is, in the anisotropic region for t < tla (spatial ex¬ 
tent: ^ < z < Id E 440 nm, where Id is the n-depletion 
region width) and in the isotropic region for t > tta {z 
> Id E 440 nm). In Fig. ll(c,d), the relative ampli¬ 
tudes of the TA modes are traced with respect to 0 in 
the isotropic and anisotropic regions, respectively. 

As can be expected intuitively in the isotropic region, 
TA mode propagation was not detected under normal in¬ 
cident probe light, as shown in Fig. 11(c); in contrast, in 
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FIG. 11: DFP oscillations with variable probe polarization at 
(a) ^ = 0° and (b) 0 = 55° under —4 V. The ^-dependent 
relative amplitudes of the TA mode, At a! Ala, in the (c) 
isotropic and (d) anisotropic region at probe incidence angles 
of 0° and 55°. 


the anisotropic region, TA modes were observed consis¬ 
tently as (j) changed, even for normal probe incidence in 
Fig. 11(d). The result reveals that the diagonal compo¬ 
nents of the dielectric tensor were perturbed not only by 
the LA modes but also by the TA modes, as indicated by 
Eq. (17). In other words, the biased region is no longer 
isotropic; consequently, both the LA and TA modes play 
roles in Maxwell’s equation for normally incident light. 
The different contributions of the modes are described in 
terms of the strain distributions along the ^ axis and the 
diagonal photoelastic coefficient for each mode. 

When the probe beam is obliquely incident on the sam¬ 
ple, the DEP interference for a TA mode is thought to 
follow the change in AR/R oc sin(0) cos(0), as denoted 
by the dashed line in Eig. 11(c) for an isotropic medium, 
because s-polarized light is flipped into p-polarization 
upon reflection from the TA wavepackets and vice versa 
for nonzero 0^. The isotropy of the unbiased region 
was confirmed by the agreement with the conventional 
sin(0) cos(0) dependence, except for a hump at 45°. The 
deviation is possibly attributed to the reduced sensitiv¬ 
ity of the LA mode due to the screening effect of the TA 
mode. The result for the anisotropic region in Eig. 11(d), 
where the light reflection is now influenced by both the 
diagonal and off-diagonal perturbations of the TA modes, 
indicates that the off-diagonal effect of the TA waves on 
the probe polarization was boosted by the amount of 
anisotropy. 

The acoustic screening effect was substantiated by the 
close relationship between the decay time tD for the LA 
mode and the relative amplitude of the TA mode sig¬ 
nal in the isotropic region, as shown in Eig. 12(a). Eor 
normal probe incidence, where the TA mode was not de- 
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FIG. 12: (a) Decay times of LA mode signals in the isotropic 
region measured with normal (squares) and 55° (triangles) 
probe incidence, (b) Decay rate of LA mode signal at ^ = 55° 
versus the relative amplitude of the TA mode. 


tected, Id was held at ^ 200 ps. On the other hand, 
the behavior oi Id was opposite to that of with 
respect to (j) ioi 0 = 55°. To further specify the rela¬ 
tionship, the decay rate l/tn for a 55° probe incidence 
is plotted as a function of in Fig. 12(b). The y 
intercept, 0.0038 ps“^, coincided with the decay rate re¬ 
sulting from the optical absorption with a skin depth of 
^2 /im. The slope, 0.025 ps“^, revealed the screening ef¬ 
fect of the TA mode on LA mode detection as AR/R = 
f dzl(ArTA)^rLA]e~^^^r)LA OC f dz[e~'^^FLA]e~^/^'>jLA, 
where the empirical decay factor 7 is 


C. External manipulation of the detection 
sensitivity in sample A 

1. Experimentally measured AC sensitivities in regions of 
different symmetry 

To gain a more quantitative insight into the relation¬ 
ship between Fi and the axial symmetry, we selectively 
measured the 5 -polarized component of the probe re¬ 
flectance as a function of 0 at —4 V, as shown in Fig. 
13. The probe polarization angle (j) was fixed at 45° to 
maximize the sensitivity of the x-polarized TA mode. As 
a result, more complicated AC propagation dynamics ap¬ 
pear in Fig. 13; the TA mode signals monotonically in¬ 
creased with 0 not only in time domain I but also in 
time domain III. The abruptly suppressed TA mode at 
tta was particularly prominent at 6>=0°. In contrast, the 
LA amplitude with respect to 0 varied in different time 
domains, slightly increasing up to tla and then rapidly 
decreasing after 2tlA‘ In an intermediate time domain, 
//, the increasing and decreasing tendencies were mixed 
for LA waves. 

The Fourier-transformed amplitudes Aj for mode in¬ 
dex i were taken in the time domains j={/,//, ///}. We 
assumed that the spatially concentrated AC strains are 
simplified to r]f=£i{z)dz • 6{z^Vit) for the ascending 

waves and rif=£i{z)dz • S(z — vA) for the descending 
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FIG. 13: Time domain signals (black line) from sample A 
for different probe incidence angles were filtered into the LA 
(blue line) and TA mode (red line) frequencies. 


waves, where £la{z) {£ta{z)) corresponds to the normal 
(shear) strain, and the region of integration, Df 
corresponds to the SDR (the i region for the LA mode 
or the n-depletion region for the TA mode). 

The ascending AC components, 77 ^, propagate toward 
the surface in time domain /, as described in Fig. 14(a); 
they then propagate back to the substrate after being 
reflected from the p-GaN/ITO interface during time do¬ 
main II in the anisotropic region, and finally enter the 
isotropic region in time domain III. On the other hand, 
the descending components, which were generated 
in the SDR, traverse the anisotropic region in time do¬ 
main I and propagate in the isotropic region toward the 
substrate when t > Ti. Therefore, A], which was exper- 


imentally extracted as shown in Eig. 
further analyzed as 

14(b-g), could be 

Al = 4Frnt + <rrdf 

{0 < t < Ti), 

Aj^ = nlrFi^'^yf + Fl^'^yf 

{Ti <t < 2rj), 

AR^ = nRrFY^yf + nR^ Ff^yf 

(2ri <t < 3Ti), 


(25) 

where Ff^'^ is the optical sensitivity of the AC 

modes for the anisotropic (isotropic) region, and r is 
the AC mode reflectivity at the p-GaN/ITO interface. 

1^1 = :^{1-e is the probe light attenua¬ 

tion factor based on the wavelength-dependent where 
j={/,//,///} = {1,2,3}. The sensitivity Ff'^^ (RV^) 
was expressed in terms of Aj {AR^) and then substituted 
into the right-hand side of the second equation in Eq. 
(25) to determine the optimal values of x = vf /{vf + yf) ^ 
which minimize the least-squares error 6{x) as 

S{x) = \Ag-Kfrai-xrii-Kj^l3i-{l-x)riif (0 < a; < 1), 

(26) 

where r]i = yf E yf. On the basis of the least-squares 
fitting by Eq. (26), x was evaluated to be ^ 0.7 (^ 0.6) 
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FIG. 14: (a) Propagation configurations of ascending ( 77 “) 
and descending (r]f) AC components in regions of different 
symmetry, (b-g) Fourier-transformed amplitudes of the DFP 
oscillations for the LA (b-d) and TA (e-g) modes in the time 
domain j={/, //, III}. 


for the LA (TA) mode, from which and F^^ could 
be experimentally extracted, as shown in Fig. 15(a-d). 


2. Electrically manipulated AC modal sensitivity 


The analytic expressions for Fi have been discussed 
separately for isotropio^ and anisotropic crystals^ in 
terms of the perturbations of the dielectric tensor due 
to strain. In isotropic materials, only the diagonal (off- 
diagonal) components in the perturbed dielectric tensor 
are induced by the LA (TA) waves. Subsequently, the 
opposite tendencies are predicted for F^^ (dashed line) 
and (dotted line) within our range of 0 values, in 

agreement with the experimental values (empty squares 
and triangles) in Fig. 15(a,b). 

On the other hand, to calculate it is necessary to 

consider the mixed nature of the dielectric tensor mod¬ 
ulations resulting from both diagonal and off-diagonal 
perturbations with empirically measured weighting fac¬ 
tors to fit the experimental results in Fig. 15 (c,d) as 


frrA ^ ^0.22 0.014\ fFn\ 
VTaJ 2.25 J \Fffy • 


(27) 


Indeed, the experimentally obtained values of Ff^'^ in 
Fig. 15(c,d) were reproduced by linear combinations 
of the diagonal perturbations in the dielectric tensor 
(oc F^X) and off-diagonal perturbations in the dielectric 
tensor (oc Most importantly, the digitized appear¬ 

ance of the typically forbidden TA mode at 6> = 0° (Fig. 
5 and 7) was theoretically verified by comparing the zero 
F^^ [Fig. 15(b)] and nonzero value of F^^ [Fig. 15(d)]. 



Incidence angle OC) 


FIG. 15: Optical sensitivity of AG modes as a function of inci¬ 
dent angle in the (a,b) isotropic and (c,d) anisotropic regions 
of sample A. 


Although further detailed investigations of the vecto¬ 
rial influence of external fields and the corresponding 
carrier-phonon interactions are necessary, it is evident 
that the mixed contributions of both diagonal and off- 
diagonal modulations in the laterally biased region re¬ 
sulted in a gradual increment of F^a increasing 

0 and allowed the TA mode detection, even at = 0^ 
(which was otherwise forbidde n^^i^^ ), demonstrating the 
electrical controllability of axial symmetry in piezoelec¬ 
tric diode structures. 


VII. SUMMARY 


In summary, we investigated the influence of the po¬ 
tential distribution on AC modal behaviors in isotropic 
structures. The crystal symmetry was broken by lateral 
electric fields in the structure with asymmetric poten¬ 
tial gradients, which enabled spatial manipulation of the 
selection rules and amplitudes of propagating acoustic 
waves. A detailed analysis of the propagation dynam¬ 
ics of the TA mode revealed abrupt suppression under¬ 
neath the laterally biased depletion region, in striking 
contrast to the electrically independent appearance of the 
LA mode. The results demonstrate the possibility of di¬ 
rect and active control of the acoustic modes in semicon¬ 
ductors simply by electrically manipulating the crystal 
symmetry. These findings could provide a novel degree 
of freedom in engineering the lattice vibrations in crys¬ 
talline heterostructures. 
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TABLE I: Matrix representations of the elastic and photoelastic tensors with different symmetries. 
Symmetry Elastic tensor Photoelastic tensor 


Hexagonal 



r^o 

^12 

^13 

0 

0 

0 ■ 

^12 

^11 

^13 

0 

0 

0 

^13 

^13 

^33 

0 

0 

0 

0 

0 

0 

044 

0 

0 

0 

0 

0 

0 

044 

0 

0 

0 

0 

0 

0 

^66. 


b°] 


Pll Pl2 Pl3 0 0 0 

Pl2 Pll P?3 0 0 0 

Pl3 Pl3 P33 0 0 0 

0 0 0 p24 0 0 

0 0 0 0 0 

0 0 0 0 0 pge 


Hexagonal with 


c=Co(l+e-,) 


• Ga 

• N 


i=ao(l-veJ 



Monoclinic with 


tan'(ej 



[n + 


n + 


^113 

^123 

^133 
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0 

0 


>fl3 

P 123 

P 133 

0 

0 

0 " 

^123 

^113 

^133 
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0 

0 


P 123 

pfl3 

P 133 
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0 

0 

^133 

^133 

^333 
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0 

0 

r-nOl A- 

P 133 

P 133 

P333 
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0 

G4^3 
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0 

^zz [P J “T 
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0 

0 

P 443 
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0 

0 

^443 

0 
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0 
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0 

0 

0 

0 

^663 J 


_ 0 
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0 

0 

0 

P663_ 
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0 

0 

0 
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0 


■ 0 
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0 

0 

P 155 

0 ■ 
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0 

0 

0 
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0 
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0 

0 

0 
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0 
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^355 
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Appendix A: Elastic and photoelastic tensors under 
different symmetries 

As introduced in Sec. IV, the elastic tensor Cim of c- 
GaN was perturbed by externally applied strains, £3 (ezz) 

and 55 ^Im ^Im ^Im 

where Ci^ is the unperturbed elastic constant, and Cimn 

is the TOE tensor—. Analogous results for the photoe¬ 


lasticity could be deduced based on the geometric sym¬ 
metry of the strained crystal 2 ii 22 as pim = + ^Pim = 

P^m Pimn^ni where p^^ is the unperturbed photoelas¬ 
tic constant, and pimn is the the tensor components for 
the nonlinear photoelasticity. In Table I, the elastic and 
photoelastic tensors under different symmetry conditions 
are comparatively summarized. 


Appendix B: derivation of formulas and 


The AC velocities for LA- and TA-mode respectively 
are obtained by solving Eq.(3) and simplified under the 
assumption that (i) > AC/m and (ii) ^ ^"^ 3 , 


vla 


ijcss + C44 + a/(C'33 - ^44)2 + 4CI5 

A 33 + ZA 7 

+ Cg^d33E, + 

Y ^33 *^44 


'^TA\\x 


^Jc 33 

+ C44 — (G33 

-C 44 ) 2 + 4 C 7|5 


V 2 po 


^ Caa 

Ch 

<^33 —<^44 





y ^^44 

+ — 

(<^355)^ fp p 2 



where Eq. and Eq. © in Sec.IV are the final form. 
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